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Abstract 

We investigate the Coulomb breakup of the two-neutron halo nucleus 6 He by extend- 
ing the framework of the finite range post form distorted wave Born approximation 
theory introduced earlier for the description of the breakup of one-neutron halo nu- 
clei. The model can use the realistic three-body wave function for the 6 He ground 
state. Calculations have been performed for energy, total and parallel momentum 
distributions of 4 He fragment. The two-body 4 He- neutron and neutron- neutron cor- 
relations have also been calculated. Our results are compared with those of an adi- 
abatic model of the Coulomb breakup reactions. The two theories lead to almost 
identical results. Comparisons are also made with the available experimental data. 
It is found that the pure Coulomb breakup can describe the bulk of the data taken 
at below grazing angles. However, a rigorous description of all the available data 
requires the consideration of the nuclear breakup effects. 
PACS numbers: 21.60.-n, 24.50.+g, 25.60.-t, 25.60.Gc 

Key words: Coulomb breakup of 6 He, two-neutron halo nucleus, post form 
DWBA with finite range effects. 



1 Introduction 



Several nuclei near the drip lines have been found to have properties which 
are strikingly different from those of their stable counterparts. These nuclei 
have a halo structure in their ground states in which loosely bound valence 
nucleon(s) has (have) a large spatial extension with respect to the respective 
core [1—5]. The two-neutron halo nuclei, of which 6 He is a prominent example, 
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have attracted a lot of attention. These are Borromean three-body systems, 
where all the two-body subsystems are unbound [3,6]. Clearly, the interaction 
of valence neutrons is vital for the stability of the two-neutron halo nuclei. 
They provide interesting examples of the three-body problem in the limit of 
very weak binding [7], and the three-body methods (within both Faddeev and 
cluster type approaches) have been extensively applied to study them [3,8-19]. 
Other approaches like non-relativistic [20] and relativistic [21,22] mean field 
models and many-body calculations [23] have also been used to investigate 
the structure of such nuclei. 

Breakup reactions, in which the valence neutrons are removed from the core in 
the interaction of halo nuclei with stable targets, provide detailed information 
on their structures [24]. The study of the breakup reaction of 6 He is interest- 
ing due to several reasons. It leads to the excitation of the continuum of this 
nucleus, which has received much attention recently due to the experimental 
discovery [25] of the soft dipole resonance (SDR), signifying the dipole oscil- 
lation between the core cluster and the valence neutrons. Although the SDR 
was predicted long ago [4,9], properties of this resonance measured in Ref. [25] 
are at variance with those predicted by the detailed three-body calculations 
[16]. The breakup reaction of this nucleus gives access to the study of the 
neutron-neutron (n-n) and neutron-alpha particle (n- 4 He) correlations, which 
are important because they provide a critical test of the dineutron model, in 
which the halo neutrons were treated as a point like "dineutron" cluster or- 
biting the core. This model has been used in earlier theoretical studies, both 
semiclassical [26] and quantum mechanical [27,28] , of the breakup of the two- 
neutron halo nuclei. Furthermore, the study of the neutron- 4 He correlation 
is expected to constrain the neutron- 4 He interaction used in the three-body 
calculations of the 6 He ground state. The study of the total and longitudinal 
momentum distributions of the charged fragment is expected to yield infor- 
mation on the matter radii of the halo. The breakup of 6 He is also important 
from the viewpoint of astrophysics as it can be used [29] to extract the rate 
of the 4 He(2n, 7) 6 He(2n, 7) 8 He reactions which could be a possible route for 
bridging the gap at the mass number, A = 5 and 8. 

Several calculations of the nuclear and Coulomb breakup [30-33] of 6 He, based 
on a model which uses a sudden approximation [34,35], have been reported 
recently. Calculations of the nuclear breakup of 6 He within the approximate 
eikonal model is reported in Ref. [36]. Due to the nature of the approxima- 
tions involved in these models, their applicability is limited to reactions at 
higher beam energies (> 200 MeV/nucleon). Moreover, they have been ap- 
plied to describe only a selected set of observables. In the distorted wave Born 
approximation (DWBA) calculations of the breakup of 6 He reported in [37], 
the interaction appearing in the kernel of the prior DWBA T-matrix has been 
evaluated in the impulse approximation which limits its applicability to higher 
beam energies. 
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In this paper, we calculate the Coulomb breakup of 6 He treating the final 
state as a four-body system. Our formalism uses the framework of the post 
form distorted wave Born approximation which is an extension of the theory 
developed earlier [24] to describe the breakup of one-neutron halo nuclei. Finite 
range of the interaction between projectile constituents is treated by making 
an effective local momentum approximation (ELMA) on the distorted wave 
of the charged core in the final channel. It was shown in Ref. [24] that this 
approximation is well satisfied for the breakup of one-neutron halo nuclei. 
Furthermore, this gave results [24] which were similar to those obtained within 
an adiabatic model (AD) [38] of Coulomb breakup reactions. This theory 
allows the use of a realistic three-body wave function for the projectile ground 
state and can be applied to breakup reactions with beam energies ranging from 
below the Coulomb barrier to higher values up to which the non-relativistic 
treatment is valid. We have used our theory to analyze almost all types of the 
breakup observables measured in the reaction of 6 He on heavy target nuclei. 
We compare our results with those of the AD and dineutron models. 

This paper is organized in the following way. Our formalism is presented in the 
next section. The structure models of 6 He are discussed in section 3. The re- 
sults of our calculations and discussions are presented in section 4. A summary 
and the conclusions of our work are given in section 5. 



2 Formalism 

We consider the breakup reaction a + 1 — > b + n\ + n 2 + 1, where the projectile 
a breaks up into fragment b (the charged core) and the two valence neutrons 
(rii and n 2 ) in the Coulomb field of the heavy target t. We consider only the 
elastic breakup mode in which the target nucleus remains in its ground state 
during the breakup process. The system of coordinates is shown in Fig. 1. The 
position vectors satisfy the relations 



— ftR, 



2m, 



(1) 



r b = r; 
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The post form T - matrix for this case is given by 
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Fig. 1. The four - body coordinate system. ni,ri2 are the valence neutrons, b and 
t represent the charged core and the target respectively. 

where c represents the cm. of the n± — n 2 system, xl-t i s the Coulomb distorted 
wave describing the motion of the charged core (b) with respect to target t 
with incoming wave boundary condition and Xc-bt describes the cm. motion 
of the valence neutrons with respect to the cm. of the b-t system. 0's are 
the internal state wave functions of the concerned particles with £'s being the 
internal coordinates. k fe and k c are the Jacobi wave vectors conjugate to r b 
and r c respectively. \I/^ + )(£ a , R, r, r^) represents the exact four-body scattering 
wave function of the projectile with an outgoing wave boundary condition. 
k a is the wave vector associated with it. V(R, r) is the sum of the two-body 
interactions between the projectile constituents, given by 

V(R, r) = V nin2 (r) + V bni (R - r/2) + V bn , 2 (K + r/2). (4) 



2. 1 Transition amplitude within the finite range distorted wave Born approx- 
imation (FRDWBA) 

In DWBA theory the breakup states are assumed to be weakly coupled and 
hence this coupling is treated only up to the first order. Thus we can write 
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R, r, ri ) = MZa, R, r)xK(ka, 



(5) 



where a (£ a ,R, r) represents the (three-body) bound state wave function of 
the projectile and Xa-t(k-a> r ») is the Coulomb distorted scattering wave de- 
scribing the relative motion of the cm. of the projectile with respect to the 
target with outgoing wave boundary condition. Substituting Eq. (5) into Eq. 
(3) and writing the latter in the integral form we get, 



rpDWBA _ J J J |^ dRdrdr . x ( : );( kb , r6 )^( 6 ) x ( : ); (kc , r j 



x 0; in2 (k, r, Un 2 )V(K, r)MU R, r)x£\(K, r t ). (6) 

We introduce the finite range effects in DWBA by an effective local momentum 
approximation by writing the distorted wave in the final channel as [24] 



X { b - > (k b ,r b ) = e* aK - R x { i ;\-k b ,r z ), (7) 

where K is an effective local momentum whose direction is taken to be the 
same as that of the asymptotic momentum k b , since the calculated cross sec- 
tions are almost independent of this quantity, as is shown in Appendix A. The 
magnitude of K is given by 



K(R d ) = ] j^(E-V(R d )). (8) 

In Eq. (8) /i b - t is the reduced mass of the b-t system, E is the energy of 
particle b relative to the target in the cm. system and V(R d ) is the Coulomb 
potential between b and the target at a distance R d . Thus, the magnitude of 
K is evaluated at a separation R d (taken to be 10 fm in all the numerical 
calculations shown in the next section), which is held fixed for all the values 
of r. We refer to Appendix A for further discussions of the validity of this 
approximation for the reactions studied in this paper. 

In Eq. (6), the Coulomb distorted waves are given by 



X £\(K,r t )=e-^T(l + t7]a y^ l F l (-tr ]a , l,t(k a r t - k a • r,)), (9) 
x£*(k 6 , r 6 ) = e *K-R e -TV2r(i + ir^e-***' 

XiFi(-i776, l,i(k b ri + k b -ri)). (10) 

For pure Coulomb breakup, the interactions between two uncharged valence 
particles as well as between their cm. and that of the b-t system are zero. 
Hence we take plane waves for the c-bt and ni-n 2 relative motions: 



5 



xi: ) fe ;(k c ,r c )=e-^-^e-^- R ! (11) 

4n 2 ( k ' r ^12)= e " lk '>n 1 n 2 ( 

Substituting Eqs. (9-12) into Eq. (6) we get, 



Cn 2 (k,r,ei 2 )=e- ik ->; tin2 (6 2 ). (12) 



rpFRDWBA = e -n( Va+Vb )/2 T ^ + -^p^ + 

x[ J dr^-^-^^i-irja, l,t(k a n - k a • n)) 
x iF^-irjb, l,i(k b ri + k 6 • r^))] 

x[ J rfrrfR^e^ (5kc - aK)R e^ k - r ^(6)0; iri2 (ei2)V(R,r) 

x a (£ a ,R,r)]. (13) 

In Eq. (13), the first integral is associated with the dynamics of the reaction 
which can be expressed analytically in terms of the Bremsstrahlung integral 
[39]. The second one (to be expressed as J^ RDWBA ^ contains the information 
about the structure of the projectile and is known as the form factor. Defining 
5k c — aK = ki, we can write 



tFRDWBA 



x R,r)]. (14) 



The function <3> a (£ a ,R, r) = V(R, r)0 o (£ a , R, r) has the following representa- 
tion in configuration space [40,41]: 

£ *tXLs(R,r)[\rt(r)®Yt(R)] L ®xs]jM. (15) 

IXLS 

In Eq. (15), the relative orbital angular momentum A between the core and 
cm. of the two neutrons is coupled with the relative orbital angular momentum 
I between the neutrons to give L, which, in turn, is coupled with the total 
intrinsic spin S of the two valence neutrons to give the resultant total angular 
momentum J of the projectile with projection M. For 6 He, J = M = 0. 



Hence 



J drdRe~^- aK ^ R e' lk - r V(R,r)(f) a ^ a ,^,r) (16) 
:47r £ (-i) l+x ® lXLS (k, fci)K(k) ® Yffyh ® Xs]oo, 



IXLS 

where 



®iXLs(k, h) = An J r 2 dr J tfdRj^kr^hR^s^r). (17) 
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The modulus square of lJ RDWBA i s given by 



jFRDW BA\2 



(4vr) 2 Yl {-) L+k ®ixLL{kM)<bvx>LL{kM)PA{™sQ) 

ll'XX'LA 

( 



x//'AA'AW(/A/'A';LA) 



I I' A 




A A' A 




(18) 



where 9 is the angle between vectors ki and k, and n = \/2n + 1. 



The dominant angular momentum configurations of 6 He contributing to Eq. 
(18) are $ 0000 (£;, k±) and $ mi (£;, hi) (with probabilities of m 0.84 and 0.11, 
respectively). Thus we have [40] 

\jFRDWBA\2 = & oooo(Kkl) + ^{kM) ~ $? m (fc,fci)P2(008e). (19) 

The differential cross section for the breakup process is given by 



d 5 a 2ir 



dE b dQ b dE ni dQ ni dfl n2 Hv a 

x \T FRDWBA \\ (20) 

where i> a is the relative velocity of the projectile in the initial channel. 

The four-body phase space factor is given by (see Appendix B) 



/ P o rp o o x _ h 9 m t m b m ni m n2 p b p ni p n2 

"^2 ~T '''712 „2 

In Eq.(21), pj's (i = t,b, 111,112) are the momenta appropriate to the frame of 
interest (i.e. laboratory or cm.) and p tot = pt + p b + p ni + p n2 - 



2.2 Transition amplitude within the adiabatic (AD) model 



The adiabatic model of Coulomb breakup assumes that the important ex- 
citations of the projectile nucleus are in the low energy continuum, i.e. the 
breakup configurations excited by the projectile-target interaction are low en- 
ergy b+n 1 +n 2 relative motion states. This allows the four-body wave function 
*i +) (£a, R, r, ri) to be written as [40] 

R, r, r t ) = MU R, r)e^ R xt\(K, r b ). (22) 
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The AD model is a non-perturbative theory and it does not make the weak cou- 
pling approximation of DWBA. However, the initial and final state Coulomb 
interactions are included to all orders in both the formalisms. It may be noted 
that while the AD model necessarily requires the valence nucleons to be charge- 
less, the FRDWBA can, in principle, be applied to cases where both the core 
and the valence particles are charged. 

The expression for the transition amplitude for the Coulomb breakup of a 
Borromean nucleus within the adiabatic model is given by [40] 



t ad = e -^ a + Vb )/i T( ^ + i7]a)T (i + i7]b) 

x [ J dr t e^' k ^> r \F l (-t Va , l,i(k a r b - k a ■ r b )) 
x iF^-irjb, l,i(k b r b + k b ■ r b ))\ 

x [|rfrrfR^e- i(kc - ak » )R e- k > 6 *(6)C in2 (62)V(R,r) 
x0 a (£ a ,R,r)]. (23) 

It is evident that except for the form factor, the transition amplitudes within 
the FRDWBA and AD formalisms are the same. The form factor in the AD 
case is given by 



x0 a (£ a ,R,r)]. (24) 

It may be noted that Eq. (23) can also be obtained within the DWBA model 
by making ELMA to the Coulomb distorted wave in the initial channel of the 
reaction, and by evaluating the local momentum at R = oo with its direction 
being the same as that of the projectile. 



3 Structure models 



The three-body wave function for the ground state (g.s.) of 6 He used in our 
calculations (which was provided to us by the authors of Ref. [17]) has been ob- 
tained [17] by solving the Faddeev equation with the hyperspherical harmonic 
expansion method which included a maximum of 20 hyperharmonics. In these 
calculations the n- A He interaction was taken from Ref. [42], and for the n-n 
interaction the Gogny-Pires-Tourreil potential [43] with spin-orbit and ten- 
sor components was used. The two-neutron separation energy and root mean 
square (rms) matter radius of 6 He predicted by this model were 0.975 MeV 
and 2.50 fm, respectively. The a-particle rms radius was taken to be 1.49 fm. 
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We would like to emphasize that in this paper our main intention is to develop 
a post form FRDWBA theory of the Coulomb breakup of the Borromean nu- 
clei where a three-body wave function for their ground state can be used. This 
enables the breakup observables to be employed as additional constraints on 
the three-body structure models of these systems. Performing structure calcu- 
lations is outside the scope of this paper. Nevertheless, it would be worthwhile 
to use the g.s. wave function of 6 He calculated by other authors (see eg. [44] 
where the calculations for the 6 He have been done within the Faddeev ap- 
proach using different sets of interactions), to see the sensitivity of our results 
on the structure models of 6 He. 

For the dineutron model calculations, a point-like 'dineutron' was assumed 
to be bound to the 4 He core in a Is state. The corresponding interaction is 
described by a Woods-Saxon potential (with radius and diffuseness parameters 
of 1.15 fm and 0.5 fm respectively) whose depth is adjusted to reproduce the 
binding energy of 0.975 MeV. The resulting 6 He wave function has the rms 
radius of 2.47 fm, with the same a-particle rms radius as that mentioned 
above. 



4 Results and discussions 

4-1 a-particle energy distribution 

In Fig. 2, we show a comparison of the calculated and measured a-particle 
energy spectra at 9 a = 5° in the elastic breakup of 6 He on a Au target at the 
beam energy of 63.2 MeV/nucleon. The data are taken from Ref. [45]. These 
cross sections have been obtained by performing the neutron angular integra- 
tions (in Eq. (20)) up to 20° and 360° in the theta and phi planes, respectively. 
The results remained almost unaffected by increasing the theta angle further. 
All the kinematically allowed values were included in the neutron energy in- 
tegration. The solid, dotted and dashed lines show the results of calculations 
performed within the FRDWBA, AD and dineutron models respectively. The 
results of the last model are plotted after multiplying the actual cross sections 
by the factor shown in the figure. It may be noted that the theoretical spec- 
trum in all the cases is shifted backwards by 10 MeV to compensate for the 
energy loss due to the target thickness [45,40]. It is seen that the FRDWBA 
results are about 20% larger than those of the AD model at the a particle 
energies near and above the peak region. However, both of them are lower 
than the experimental data. The grazing angle for the reaction studied in this 
figure is about 3.5°. Therefore, nuclear breakup effects are expected to make 
significant contribution to the cross sections at the angle for which these data 
are taken. This could explain the pure Coulomb breakup results being smaller 
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Fig. 2. Calculated energy distributions of a-particle in the elastic breakup of 6 He on 
a Au target at 63.2 MeV/nucleon at 9 a = 5°. The solid, dotted and dashed curves 
represent the results of the FRDWBA, AD and dineutron models respectively. The 
results of the di-neutron model have been plotted after multiplying the actual cross 
sections by 0.33. The data have been taken from [45]. 



than the experimental cross sections, particularly in the peak region. 

It may be noted that the 'dineutron' model overestimates the measured cross 
sections by a factor of about 3 (the dashed line). This result can be understood 
qualitatively by recalling that the magnitudes of the calculated cross sections 
are dependent on the mean square distance between the a-particle and the 
center of mass of 6 He [46]. In case of the three-body wave function, this dis- 
tance is 1.428 fm 2 , whereas for the 'dineutron' model, it is 2.32 fm 2 . The larger 
cross section obtained in the latter model, therefore, may be attributed to this 
difference. 

Measurements performed below the grazing angle are expected to be less af- 
fected by the nuclear breakup. One could then hope to test the significance 
of the small but noticeable difference seen in the pure Coulomb breakup cross 
sections calculated within the FRDWBA and AD models. However, it can 
already be seen that the peak position in the cross section calculated within 
the former is somewhat shifted to smaller energies as compared to that of the 
latter and is in a slightly better agreement with the peak position of the data. 
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Fig. 3. Calculated Coulomb part of total two- neutron removal cross sections for 
breakup of 6 He on a Pb target at three different beam energies of 30, 40 and 50 
MeV/nucleon. The open triangles, open circles and open squares show the results of 
the AD, FRDWBA and dineutron models respectively. The data point is the average 
of the experimental total two-neutron removal cross sections measured between the 
beam energies of 28 MeV/nucleon to 52 MeV/nucleon ([47]). 

4-2 Coulomb part of total two-neutron removal cross sections 



In Fig. 3, we show the results of calculations for the Coulomb part of the total 
two-neutron removal cross section at three different beam energies of 30, 40 
and 50 MeV/nucleon for the breakup reaction of 6 He on a 208 Pb target. The 
open circles, open triangles and open squares show the results of the FRD- 
WBA, AD and dineutron models respectively. With increasing beam energy 
the cross sections of first two models decrease, which is consistent with the 
earlier observations (see e.g. [27]). We note that here also the results of these 
two models differ from each other by 20-25%, but both of them are smaller 
than the rough estimates for the pure Coulomb breakup cross sections given 
in Ref. [47]. However, such estimates are unlikely to be accurate for nuclei 
having a halo structure as discussed in [48,49]. 

That the nuclear breakup contributions are substantial in the case of breakup 
of 6 He even on a heavy target, is evident from this figure, where we also show 
(solid circle) the experimental value of the measured total two-neutron removal 
cross section (crf. 3 ^), which is the average of the measurements performed at 
beam energies between 28 to 52 MeV/nucleon. Our theoretical values corre- 
spond only to the Coulomb part of the elastic breakup cross section, while 
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<j e ^2n has contributions from the Coulomb, nuclear as well as their interference 
terms. Moreover, the nuclear breakup consists also of the inelastic breakup 
mode [50] (which is called as the stripping process in Ref. [36]). Therefore, 
to understand the data of Ref. [47], it is essential to calculate the nuclear 
breakup cross sections on the same footing as the Coulomb one. Similarly, 
in the breakup of 6 He on the Pb target at the higher beam energy of 240 
MeV/nucleon, the calculated pure Coulomb two-neutron removal cross sec- 
tion is 62.5 mb, which is quite small as compared to the corresponding value 
of a e *2 n [51]. This difference too can be understood from the fact that the cross 
sections corresponding to the nuclear elastic and inelastic breakup modes are 
quite large [33,37]. 

The pure Coulomb two-neutron removal cross sections, calculated within the 
dineutron model (open squares), are even larger than crf^, and have a wrong 
type of the beam energy dependence. Thus, the dineutron model is not suitable 
for describing the data for the breakup of 6 He. 

4-3 a-particle momentum distributions 

8000 i i i 




1100 1300 1500 1100 1300 1500 

p a (MeV/c) 

Fig. 4. Comparison of the calculated and measured total momentum distributions 
of the alpha particle in the breakup of 6 He on a Au target at 64 MeV/nucleon at 
four forward angles. The results of the FRDWBA and AD models are shown by the 
solid and dotted lines, respectively. The data are taken from [52]. 

In Fig. 4, we compare the calculated and measured [52] total momentum 
distributions of the a-particle emitted in the breakup of 6 He on the Au target 
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at the beam energy of 64 MeV/nucleon. The ct-particle was detected at four 
average angles of 0.65°, 1.67°, 3.02° and 4.55° [52]. In this experiment, the 
coincident 7-ray technique was used to identify events in which the target 
nucleus was simultaneously excited. The elastic breakup events were tagged 
with 7-ray multiplicity M 7 = 0. Since the data are given in arbitrary units, we 
have normalized our calculations to the peak of the experimental distributions 
in all four cases. It is clear that shapes of the momentum distributions are well 
reproduced by our calculations. The full width at half maximum (FWHM) of 
the momentum distributions changes systematically from 76 MeV/c to 98 
MeV/c with increasing angle. This could indicate that the effects related to 
the reaction dynamics (e.g. distortion of the wave functions describing the 
fragment-target relative motion) lead to the widening of the total momentum 
distributions as a function of the detection angle. 

It should be remarked here that the calculations [53], where the total mo- 
mentum distributions are obtained from the momentum space wave functions 
of the ground state of the projectile disregarding the reaction dynamics al- 
together, are unlikely to describe this feature of the experimental data. At 
the same time, if the detection angle increases further the nuclear breakup 
effects will become more important and pure Coulomb breakup may also be- 
come inadequate to describe the FWHM of the experimental total momentum 
distributions. Therefore, it would be interesting to span a larger range of de- 
tection angles in future experiments. 

30 
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Fig. 5. Parallel momentum distribution of the alpha particle emitted in the Coulomb 
breakup of 6 He on U and Pb targets at the beam energy of 24 MeV/nucleon. The 
FRDWBA and AD model results are shown by the solid and dotted lines respec- 
tively. The data are taken from [54] . 
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In Fig. 5, we present the comparison of our calculations with the data for the 
parallel momentum (p 2 ) distribution of the a-particle emitted in the breakup 
reaction of 6 He on U and Pb targets at the beam energy of 24 MeV/nucleon. 
These measurements have recently been performed at the Michigan State Uni- 
versity [54]. The calculated distributions have been obtained by integrating 
Eq. (20) over the neutron angles up to 5.7° (which is the minimum opening 
angle of the neutron detector in this experiment [54]). The calculated results 
are normalized to the peak of the data. We note that although both FRD- 
WBA and AD models are able to reproduce the FWHM of the experimental 
p z distributions (~ 76 MeV/c) reasonably well, the former is in slightly better 
agreement with the data. 

4-4 Neutron- neutron and core-neutron relative energy spectra 




Fig. 6. Calculated neutron-neutron relative energy spectra in the Coulomb breakup 
of 6 He on the Pb target at the beam energy of 240 MeV/nucleon. The FRDWBA 
and adiabatic model results are shown by the solid and dotted lines respectively. 
The data are taken from [51]. 

In Fig. 6, we show the comparison of our calculations with the data [51] for 
the neutron-neutron relative energy (E nn ) spectrum obtained in the breakup 
of 6 He on a Pb target at the beam energy of 240 MeV/nucleon. The calcu- 
lated spectrum has been obtained by transforming the cross section shown in 
Eq. (20) to the frame of relative and cm. motion of the fragments (see Ap- 
pendix B), and integrating the resultant over the unobserved variables. The 
integrations over the theta and phi angles of n\-U2 and a-(ni-n 2 ) systems were 
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carried out in the range of 0°-180° and 0°-360° respectively, while those over 
the energy E Q _( ni _ n2 ) are done in the range of 0-10 MeV. However, the in- 
tegration over the theta angle of the 6 He-target system (6 at in Eq. (B.12)) is 
carried out up to the grazing angle (1.1°) only, as we are interested only in 
the pure Coulomb contributions to the cross sections. The solid (dotted) line 
shows the corresponding FRDWBA (AD) model results. Both calculations are 
normalized to the peak of the experimental data which is at E nn ~ 0.45 MeV. 
Both the models lead to almost similar results which incidently are in good 
agreement with the shape of the experimental data. Of course, inclusion of 
the nuclear and Coulomb-nuclear interference terms can affect both the ab- 
solute magnitude as well as shape of this distribution. Therefore, it would be 
interesting to have the data with absolute values on the one hand, and the 
calculations including both Coulomb and nuclear breakup effects on the other. 

The shape of the neutron-neutron energy spectrum reported in Ref. [55] is 
not in agreement with that of the experimental data. This is due to the fact 
that these authors have considered a pure sudden approximation, and have 
ignored the final state interaction effects. This underlines the importance of 
the reaction dynamics in the breakup reactions even for beam energies around 
240 MeV. 




E an (MeV) 

Fig. 7. Calculated n-a relative energy spectra in the Coulomb breakup of 6 He on 
the Pb target at the beam energy of 240 MeV/nucleon. The FRDWBA and AD 
model results are shown by the solid and dotted lines respectively. The data are 
taken from [51]. 

In Fig. 7, we have compared the calculated n- 4 He relative energy (E an ) spec- 
trum with the corresponding experimental data [51]. The peaks of the theo- 
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retical calculations are normalized to the maximum of the experimental data. 
The theoretical cross sections have been obtained by integrating Eq. (B.14) 
of Appendix B over the unobserved variables. The range of the angular and 
the energy integrations of various variables were the same as that described 
above in connection with Fig. 6. We see that except for the peak position, the 
theoretical calculations reproduce the data reasonably well. The difference in 
the peak positions of the experimental and theoretical spectra may perhaps 
be attributed to the fact that the position of the 5 He resonance is not properly 
reproduced within the three-body model used in getting the ground state wave 
function of 6 He employed in our calculations. We plan to check this point in 
a future study by using other 6 He wave functions which may be better in this 
respect. 



4-5 Excitation energy spectra 




Fig. 8. Calculated excitation energy spectra in the Coulomb breakup of He on the 
Pb target at the beam energy of 240 MeV/nucleon. The FRDWBA and AD model 
results are shown by he solid and dotted lines respectively. The data are taken from 
[51]. The dashed line shows the result of the dineutron model. 

In Fig. 8, we show the pure Coulomb contribution to the excitation energy 
(E ex ) spectrum in the breakup of 6 He on a Pb target at the beam energy of 240 
MeV/nucleon [51]. Excitation energy is defined as the sum of the a-(ni-n 2 ) 
and ni~n 2 relative energies and the Q-value of the reaction. The calculated 
cross sections are normalized to the peak of the data as the elastic Coulomb 
breakup contribution constitutes only a part of the excitation energy spectrum 
[37]. We note that both FRDWBA and AD models give identical results and 
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are able to reproduce the shape of the spectrum well. We have also shown 
in this figure the results obtained with the dineutron model, which is clearly 
inadequate to describe the data. It should, however, be noted that the cal- 
culated cross sections have not been convoluted with the detector response 
matrix as suggested in [51]. As is shown in [51], this folding reduces the cross 
sections by about 20 % near the peak region, and by less that 5% beyond it. 
Therefore, the agreement with the data seen in Fig. 8 should be seen in this 
context. 



5 Summary and Conclusions 

In this paper we have calculated the elastic Coulomb breakup of the Bor- 
romean nucleus 6 He within a four-body post form distorted wave Born ap- 
proximation approach. The finite range effects have been included by using an 
effective local momentum approximation on the Coulomb distorted wave func- 
tion of the charged heavy core in the outgoing channel. This method allows 
the use of realistic correlated three-body wave functions for the projectile. 
Thus the special features associated with the structure of the two-neutron 
halo systems can be taken into account in the reaction calculations. In our 
method, the breakup amplitude is expressed as a product of the factors de- 
scribing separately the projectile structure and the dynamics of the reaction. 
This feature of our theory is common with that of the adiabatic model of the 
Coulomb breakup reactions. Although unlike DWBA, this model does not use 
the weak coupling approximation to describe the center of mass motion of 
the fragments with respect to the target, yet it leads to results which agree 
with those of the FRDWBA within approximately 20%. The existing data are 
unable to highlight the significance of this small difference between the two 
calculations. Measurements done below the grazing angle may be helpful in 
this respect. 

We have presented calculations for the energy, total and parallel momentum 
distributions of the a-particle as well as for the neutron- neutron and a-neutron 
relative energy spectra. In all the cases, our calculations are able to reproduce 
the shapes of the experimental spectra reasonably well. However, the abso- 
lute magnitudes of the data (wherever available) were underpredicted by both 
finite range DWBA and adiabatic model of the Coulomb breakup. This sug- 
gests that nuclear breakup contributions are substantial in the kinematical 
regimes covered by these data. Furthermore, although the pure Coulomb cal- 
culations are able to reproduce the shapes of the experimental data on the 
fragment-fragment correlations, inclusion of nuclear breakup effects are neces- 
sary for a more comprehensive description of the data. Our calculations rule 
out the point-like dineutron-core type of models for describing the breakup 
of Borromean systems. The two-body fragment-fragment correlation data are 
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beyond the scope of such models. 
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A Validity of the effective local momentum approximation 
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Fig. A.l. Variation of rj(r) (upper half) and K(r) (lower half) with r for the reaction 
6 He + 197 Au -► 4 He + n + n + 197 Au at the beam energies of 240 MeV/nucleon 
(left side) and 63.2 MeV/nucleon (right side). 



The derivation of the finite range amplitude [Eq. (13)] uses an effective local 
momentum approximation (ELMA). This allows us to write the distorted 
wave function of argument product of the plane wave of argument R 

and distorted wave of argument r,. This simplifies the evaluation of the six- 
dimensional integrals by reducing them as a product of two three-dimensional 
integrals of independent arguments. 
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As discussed in [56], a condition for the validity of this approximation is that 
the quantity 

v ( r ) = -KWldKirf/drl' 1 (A.l) 

evaluated at a representative distance Rd should be much larger than the rms 
radius of the projectile. To check this, we show in Fig. A.l, r](r) (upper half) 
and K(r) (the magnitude of the local momentum) (lower half) as a function 
of r, for the breakup reactions 6 He + 197 Au — > 4 He + n + n + 197 Au at the 
beam energies of 240 MeV/nucleon (left side) and 63.2 MeV/nucleon (right 
side). We see that for r > 8 fm, i](r) is several orders of magnitude larger than 
rms radii of the halo in both the cases. Therefore, the above condition is well 
satisfied. 

/From the lower half of Fig. A.l, we note that the value of K(r) remains 
constant for distances larger than 10 fm for both the beam energies. Due to the 
peripheral nature of the breakup reaction, this region contributes maximum 
to the cross section. Therefore, our choice of a constant magnitude for the 
local momentum evaluated at 10 fm is well justified and is independent of the 
beam energy. In fact, we noted that as Rd is increased from 5 to 10 fm the 
calculated cross sections vary by at the most 10 %, and with a further increase 
the variation is less than 1 %, in all the cases considered in this paper. 

In the application of the ELMA in the description of the heavy ion induced 
transfer reactions [57] and the breakup reaction of the one-neutron halo nu- 
clei [24], the calculated cross sections were found to be almost unaffected by 
the choice of the direction of the local momentum. To study the sensitivity 
of our FRDWBA results on the direction of K for reactions under study in 
this paper, we have performed calculations for three cases where we take the 
angles of the local momentum {d\) parallel to those of k&, (d 2 ) parallel to the 
direction corresponding to the half of the angles of and (c^) parallel to the 
beam direction (zero angles). In the upper part of Fig. A. 2, we show the energy 
distribution of the alpha particle for the same reaction as described in Fig. 2. 
The results obtained with cases d±, d 2 and d 3 are shown by solid, dashed and 
dotted lines respectively. We note that energy distributions calculated with 
these choices are almost identical. The results for the neutron-neutron energy 
correlations for the same reaction as in Fig. 6, are shown in the lower part 
of this figure for the choices (di), (d 2 ) and (<i 3 ) for the direction of the local 
momentum. In this case too we observe almost no variation in the calculated 
cross sections. Similarly, we have noted no dependence of the calculated widths 
of the parallel momentum distributions of heavy fragments on the direction of 
K, in all the reactions investigated in this paper. Therefore, the dependence of 
various cross sections for the reactions studied in this paper, on the direction 
of the local momentum is either very minor or almost negligible. The mea- 
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E nn (MeV) 



Fig. A. 2. The energy distributions of the alpha particle in the reaction He + Au 
— ► 4 He + n + n + 197 Au at the beam energy of 63.2 MeV/nucleon (upper part), 
and the neutron-neutron energy correlation in the reaction 6 He + 208 Pb — > 4 He + 
n + n + 208 Pb at the beam energy of 240 MeV/nucleon (lower part), for choices 
d\ (solid line), d,2 (dashed line) and d% (dotted line) for the direction of the local 
momentum, as discussed in the text. 

surements done so far are not able to distinguish the small differences that we 
see here in some cases. Hence, we have performed all our calculations in this 
paper by using the direction of the momentum K to be the same as that of 



B The four body phase space factor for a given Q value 

The phase space in which we fix the total energy E tot in the final state at the 
value E is given by 

da = h~ 9 5(E tot - E)dp b dp ni dp n2 (B.l) 
The total energy is given by 

^ = E a + Q (B.2) 
= E t + Eb + E ni + E n2 (B.3) 
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■.E b + E ni + 7 ^ + ^- (B.4) 

2m„ 2 2m t 

= £ 6 + £ m + + ^-[(Ptot -Pb~ PnJ 2 + p\ 2 (B.5) 

2 Pn 2 -(Ptot - P& - P 



where i? a is the projectile energy and Q is the Q- value of the reaction (<0). 
Ptot is the total momentum of the projectile in the incident channel, which by 
momentum conservation, is given by 

Ptot = Pt + P6 + Pni+Pna- ( B - 6 ) 

Using dp = m.p.dE.dQ, we have 



da da ^ 

dE b dVl b dE ni dVL ni dE n2 dVt n2 dp b dp ni dp n2 

= h~ 9 5(E tot - E)m b p b m ni p ni m n2 p n2 (B.8) 

Therefore 



p(E b , Q b , E ni , Vt ni , £l n2 ) 



da 

dE b dVl b dE ni dVl ni dVl n2 
da 



-dE, 



dE b dVt b dE ni dVl ni dE n2 dVl n 
h' 9 m b m ni m n2 J p b p ni p n2 5 (E tot - E)dE, 



= h 9 m b m ni m n2 



PbP ni Pn 2 



dE tot /dE,, 



Tl'2 



Etot = E 



(B.9) 



Now from Eq.(B.5) 



dE tot dEtot dp n2 1 Pn 2 -(Ptot-Pb 

~ry^i = "5 -JW~ = m "2 + m t- m n 2 o 

dE n2 dp n2 dE n2 m t p z n2 
Therefore 



(t? o p o o ^ h 9m t m bm ni m n2 p b p ni p n2 

p{E b ,U b ,E ni ,U ni ,U n2 ) = ■ Pn9 .( Ptot - Pb - Pn J - ( B - n ) 

m ri2 +m t - m ri2 2 , lJ - 

fn 2 



It can also be shown that 
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d 5 a 



dE rilTl2 dVl nin2 dE b —(j lin2 ^dVl b ^(, niri2 ^dVt a t 



x ^1 , (B.12) 

dE b dVt b dE ni dVl ril dVt n2 



where 



j m n 1 n2Pn 1 n2 m b—(nin2)Pa-{n\n2) m atPat 

and 



J- 



2 



dE bni dVL bni dE n2 _ ( 6ni ) df2„ 2 _ ( 6ni ) c/S7 at 

x 

dE b d£l b dE ni dVt ni dVl n , 2 



x Jc . ^ ■ (B.14) 



where 



^ _ fn bni p bni m n2 _ ( 6ni )p TO2 _ ( 6ni ) m at p at 

■> E ni , fi ni , ) 
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